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ABSTRACT 

We present VLA detections of radio emission in four four-image gravitational lens 
systems with quasar sources: HS 0810+2554, RX J0911+0511, HE 0435—1223 and 
SDSS J0924+0219, and e-MERLIN observations of two of the systems. The first 
three are detected at a high level of significance, and SDSS J0924+0219 is detected. 
HS 0810+2554 is resolved, allowing us for the first time to achieve 10-mas resolution of 
the source frame in the structure of a radio quiet quasar. The others are unresolved or 
marginally resolved. All four objects are among the faintest radio sources yet detected, 
with intrinsic flux densities in the range 1-5/iJy; such radio objects, if unlensed, will 
only be observable routinely with the Square Kilometre Array. The observations of 
HS 0810+2554, which is also detected with e-MERLIN, strongly suggest the presence 
of a mini-AGN, with a radio core and milliarcsecond scale jet. The flux densities of 
the lensed images in all but HE 0435—1223 are consistent with smooth galaxy lens 
models without the requirement for smaller-scale substructure in the model, although 
some interesting anomalies are seen between optical and radio flux densities. These 
are probably due to microlensing effects in the optical. 
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1 INTRODUCTION 

Strong gravitational lens systems, in which background 
sources are multiply imaged by foreground galaxies, are 
important for two main reasons. Eirst, the lensing effect 
magnifies the background source; although the surface 
brightness is conserved, the area increases and we can 
observe background sources with typically 5-10 times 
better signal-to-noise than without the lensing. Second, 
the lensing effect allows us to probe the mass distri¬ 
bution of the lensing galaxy, on scales from the over¬ 
all mass profile down to the scales of individual stars. 
General reviews of strong lensing and its applications 
relevant to this work are given by Kochanek (2004), 
Gourbin, Saha & Schechter (2002), Zackrisson & Riehm 
(2010), Bartelmann (2010) and Jackson (2013). 

Lens systems in which the background source is 
a quasar were the first class of systems to be discov¬ 


ered, mostly in radio surveys (Walsh, Garswell & Wey- 
mann 1979; Hewitt et al. 1988; Browne et al. 2003) but 
later in optical surveys (e.g. Wisotzki et al. 1993; In- 
ada et al. 2003). Lens systems with optical, or “radio¬ 
quiet” , quasars as the source now dominate the sample 
of strongly lensed quasars. The radio-selected sample is 
mostly composed of the 22 lenses from the Gosmic Lens 
All-Sky survey (GLASS; Myers et al. 2003; Browne et 
al. 2003) together with smaller samples from the MIT- 
Greenbank (MG) and southern surveys (e.g. Hewitt et 
al. 1992; Winn et al. 2002). This sample has not ex¬ 
panded significantly in the last 15 years, because of the 
difficulty - in the era before the Square Kilometre Array 
- of conducting more sensitive wide-fleld radio surveys 
at the necessary sub-arcsecond resolution for lens dis¬ 
covery. 

In order to increase the size of the sample of quasar 
lenses observed at radio wavelengths, we can use the 
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fact that all quasars are likely to have radio emission at 
some level. For example, White et al. (2007) performed 
a stacking analysis at the positions of radio-“quiet” 
quasars from the FIRST 20-cm radio survey (Becker 
et al. 1995) and found that quasars not detected at the 
1-mJy level are likely, on average, to have radio flux 
densities which fall with decreasing optical brightness 
to ^ 100/iJy at / = 20 — 21. This is within reach of the 
new generation of telescopes such as the Karl G. Jan¬ 
sky Very Large Array (VLA) and the extended Multi- 
Element Remote Linked Interferometer (e-MERLIN). 
We began a programme (Wucknitz & Volino 2008, Jack- 
son 2011) to investigate four-image quasar lens sys¬ 
tems without current radio detections, and achieved a 
successful detection, at around the 20-30/iJy level, of 
lensed images of the background quasar in the cluster 
lens system SDSS J1004+4112, as well as detections 
of lensed images at a brighter level in the lens system 
RXJ 1131-1231. 

There are three main motivations for radio observa¬ 
tions of the radio emission from radio- “quiet’Q quasar 
lens systems. The first is that the emission mechanisms 
of radio-quiet quasars are not well understood, and in 
particular it is not clear how far the mechanisms which 
power radio-loud quasars - accretion and the forma¬ 
tion of a powerful jet close to a black hole - also ap¬ 
ply to radio-quiet quasars. In sources of intermediate 
flux density, jets characteristic of ordinary active nu¬ 
clei are seen (Blundell & Beasley 1998, Leipski et al. 
2006) and the radio sources appear to be variable (Bar- 
vainis et al. 2005). Both of these observations suggest 
that the AGN paradigm applies to these objects. In 
fainter cases (<1 mJy at 1.4 GHz at redshifts of a few 
tenths) however, Gondon et al. (2013) argue that star 
formation is the primary mechanism for the emission. 
This is inferred from an analysis showing that the ra¬ 
dio luminosity function of QSOs turns up sharply below 
Ti. 4 ghz < 10^^ WHz“^ and suggesting that a second 
population emerges at these luminosities. It is also pos¬ 
sible that a different emission mechanism is at work in 
the core of the sources. Blundell & Kuncic (2007) sug¬ 
gest the presence of optically-thin bremsstrahlung emis¬ 
sion (although see Steenbrugge et al. 2011 for evidence 
against this view) and Laor & Behar (2008) propose the 
possibility of emission from magnetically-heated coro- 
nae, rather than a classical AGN source. It is important 
to achieve high-resolution radio imaging of these sources 
in order to separate the possibilities. This is typically 
very difficult to achieve with current long-baseline in¬ 
terferometers in very faint objects. However, the use of 
leasing magnification provides a way to detect otherwise 
unobservable objects. 

The second motivation is that radio and optical ob¬ 
servations are subject to different propagation effects. 
The main such effect in the optical is microlensing due 
to stars in the lens galaxy, which produces measurable 
flux changes because the size of the optical source is 


^ In the rest of this article, we drop the inverted commas, but 
it should be understood that radio-“quiet” quasars are not radio 
silent. 


very small. Repeated optical monitoring can reveal the 
flux density variations associated with the movement 
of the source with respect to the microlensing caus¬ 
tic patterns (Irwin et al. 1989, Wisotzki et al. 1993; 
Burud et al. 2002; Poindexter, Morgan & Kochanek 
2008; Blackburne et al. 2011; Munoz et al. 2011). At 
radio wavelengths, the source is larger, and therefore 
much less susceptible to microlensing; hence, compari¬ 
son between the two wavebands can reveal the effects 
of microlensing directly. In the radio, the only signifi¬ 
cant propagation effect is scattering by ionized media 
(Koopmans et al. 2003This seems to be particularly 
noticeable in a few cases, such as CLASS B0128+437 
(Phillips et al. 2000, Biggs et al. 2004) but should in 
principle decrease strongly at higher radio frequencies. 
A corresponding problem at optical wavelengths is ex¬ 
tinction by dust in the leasing galaxy, which can be 
used to learn about the properties of the dust if the 
intrinsic fluxes are known (e.g. Jackson, Xanthopoulos 
& Browne 2000, Eliasdottir et al. 2006, Ostman et al. 
2008). Radio wavelengths therefore provide an impor¬ 
tant input to any programme which aims to disentangle 
the effects of substructure in the lens galaxy from those 
of microlensing and extinction. 

The third motivation for radio observations of four- 
image gravitational lens systems is their potential to 
probe substructure in the lens galaxies. In principle, 
quasar lens systems are useful for probing small-scale 
structure within the lens potential, down to 10 ^Mq and 
below (Mao & Schneider 1998), because the flux den¬ 
sity of the lensed image can be perturbed by small-scale 
mass structures close to the corresponding ray path. 
Such sub-galactic scale structures are important pre¬ 
dictions of Cold Dark Matter (CDM) models. Initially 
they were thought not to be present in required quanti¬ 
ties in the Milky Way (Moore et al. 1999, Klypin et al. 

1999) . The situation is now less clear, as a population 
of faint Milky Way satellites have since been discovered 
(Belokurov et al. 2006; Zucker et al. 2006; Koposov et 
al. 2015). The Milky Way halo mass is a critical variable 
(Wang et al. 2012; Kafle et al. 2014) as the predicted 
halo population is sensitive to it. In lens systems, the 
flux density of the lensed images is particularly sensitive 
to small structures, because it depends on the second 
derivative of the leasing potential, as opposed to the 
image positions which depend on the first derivative. 
The usual evidence for a detection of substructure is 
therefore a set of image flux ratios which cannot be fit 
by smooth models. More particularly, four-image lenses 
in cusp configurations (where the source is close to the 
cusp of the astroid caustic) and fold configurations (re¬ 
sulting from the source being close to the caustic fold) 
give clear theoretical expectations for image flux ratios 

^ In principle, the size of a compact, synchrotron self-absorbed 
radio source decreases as the square root of the brightness, but 
this is unlikely to result in microlensing until the source is fainter 
than I/rJy. Claims exist for radio microlensing (Koopmans et al. 

2000) which could also be explained by other propagation effects. 
In principle, free-free absorption is also possible (Mittal et al. 
2007) although the electron columns are likely to be too small in 
all but exceptional cases. 
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which must be obeyed by any smooth model. In cusp 
lenses, there are three close images and the middle im¬ 
age is expected to have the brightness of the sum of 
the outer two (Schneider & Weiss 1992); in fold lenses, 
the two close images are expected to have the same flux 
(Keeton, Gaudi & Fetters 2003; Congdon, Keeton & 
Nordgren 2008). Because of optical microlensing, radio 
lens systems have traditionally been used for this work 
(Mao & Schneider 1998, Fassnacht et al. 1999, Met¬ 
calf & Zhao 2002, Metcalf 2005, Kratzer et al. 2011). 
The statistics of such objects were analysed by Dalai & 
Kochanek (2002) who found a fraction of between 0.6% 
and 7% of mass in 10^ — 10^ Mq substructures at the 
radius probed by the lensing. More recent theoretical 
work has used more realistic treatment of lens galax¬ 
ies via the use of numerical simulations (Bradac et al. 
2004, Amara 2006, Maccib et al. 2006, Xu et al. 2009). 
The conclusions are generally that there is, if anything, 
an excess of substructure over what is predicted to ex¬ 
ist in CDM (though see Metcalf & Amara 2012, Xu et 
al. 2015). At the same time, analyses of individual lens 
systems have yielded constraints on substructures at the 
^ IO^Mq level (e.g. Fadely & Keeton 2012). The major 
problem in using quasar lenses to constrain substruc¬ 
ture has been the lack of large enough samples of radio- 
loud four-image lenses; the seven studied by Dalai & 
Kochanek in 2002 have formed the sample for most sub¬ 
sequent investigations. Possible solutions to this prob¬ 
lem include the use of mid-infrared fluxes (Chiba et al. 
2005; Fadely & Keeton 2011) assuming that the mid-IR 
source is large enough not to be subject to microlensing 
(but see Sluse et al. 2013). An alternative approach is 
to use the narrow line region of quasars (Moustakas & 
Metcalf 2003, Sugai et al. 2007, Nierenberg et al. 2014) 
which should also be large enough to be unaffected by 
microlensing, or submillimetre observations in the case 
of new lenses from Herschel and the SPT (Hezaveh et 
al. 2013). A further alternative is to perform lens re¬ 
construction of systems with extended sources (War¬ 
ren & Dye 2003, Dye & Warren 2005, Koopmans 2005, 
Vegetti et al. 2009, Vegetti et al. 2012). Cases of sub¬ 
structure detections have already been reported from 
these works, although the sensitivity is mainly towards 
the higher-mass end of the substructure mass function; 
quasar lenses are thus likely to be usefully complemen¬ 
tary to this method. 


This work presents a continuation of a programme 
to detect and image faint radio sources in gravitational 
lens systems. Its aim is to increase the number of four- 
image lenses with detected radio fluxes, both to increase 
the sample sizes of quasar lens systems suitable for the 
investigation of sub-galactic scale substructures in the 
lens, and to begin the study of the very faint radio 
sources which are imaged by the foreground lens galax¬ 
ies. Where necessary we use a standard flat cosmology 
with Qrn = 0-27 and Hq = 68kms“^Mpc“^. 


2 SAMPLE AND OBSERVATIONS 
2.1 Sample selection 

Our target sample includes all known gravitational lens 
systems with four lensed images, no detected radio emis¬ 
sion above the ^1 mJy level reached in large-scale sky 
surveys such as the FIRST and NVSS 1.4-GHz surveys 
(Becker et al. 1995, Condon et al. 1998), and with dec¬ 
lination > —20° for accessibility to the VLA and e- 
MERLIN radio arrays. There are 13 of these in cur¬ 
rent compilations such as the CASTLES (Kochanek et 
al. 1998) and Masterlens (Moustakas et al. 2012) lists, 
which represents a potential factor of 3 improvement 
in statistics if radio flux densities can be measured for 
all of them. One of these objects, SDSS J1004+4112, 
was already detected by Jackson (2011) using the VLA 
in the lower resolution C-configuration. This is a wide- 
separation object (Inada et al. 2003a) resulting from 
the lensing action of a cluster, whose mass distribution 
is correspondingly more difficult to model. Most such 
objects, however, are lensed by individual galaxies; we 
have in this preliminary phase used the VLA (resolution 
^073 at 5 GHz) to investigate the wider-separation ob¬ 
jects within this sample. A further object, RXJ 1131- 
1231, was previously detected by Wucknitz & Volino 
(2007) in archival VLA data, and subsequently with the 
VLA and MERLIN (although not with VLBI). Table 
shows a summary of the lens systems observed (inclucF 
ing, for completeness, SDSS J1004+4112) together with 
the source and lens redshifts and other observational in¬ 
formation. 


2.2 Observations and data reduction 

2.2.1 VLA observations 

Objects were observed with a total bandwidth of 2 GHz 
in 16 lEs of 128 MHz over the frequency range 4488- 
6512 MHz (C-band). The exception was HS 0810+2554 
which was observed at X-band, with a similar spectral 
arrangement over the frequency range 7988-10036 MHz, 
in order to achieve the resolution needed for the smaller 
spatial scale of this source. Integration times were gen¬ 
erally a few hours (Table [B and observations were taken 
at various times during the autumn of 2012. Although 
the individual 3 or 1.5-hour observations were taken 
at different times, we do not detect variability in any 
case between individual epochs. All new observations 
were taken in A-configuration, which has a maximum 
baseline of 36 km and a theoretical resolution of 0735 
at 5 GHz and 0722 at 8.4 GHz. Data were taken in 
spectral-line mode, with 2-MHz channels, although this 
was reduced in subsequent processing as only a small 
area of sky was required. 

Nearby phase calibrators were observed at regular 
intervals, with a pattern of 1 minute on the calibrator 
and 5 minutes on source, in order to calibrate the in¬ 
strumental and atmospheric phases. Sources of known 
flux density, either 3G138 or 3G286, were observed in 
order to fix the flux density to the Baars et al. (1977) 
scale. 

Data analysis was performed using the NRAO AlPS 
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Object 

-^-lens 

-S^source 

Separation 
/ arcsec 

References 

VLA Obs. date 
(2012) 

Frequency 

/GHz 

Exposure 

time/h 

noise 

/fjJy/h 

HE 0435-1223 

0.46 

1.69 

2.5 

W02,M05,006 

Oct 26, Nov 9 

5 

6 

3 

HS 0810+2554 

? 

1.50 

0.9 

R02 

Oct 22, Nov 8,24 

8.4 

7.5 

3 

RX J0911+0551 

0.77 

2.80 

3.2 

B97,B98,K00 

Oct 31. Nov 6,24 

5 

7.5 

3 

SDSS J0924+0219 

0.39 

1.52 

1.8 

I03B,E06,006 

Nov 5 

5 

3 

3 

SDSS J1004+4112 

0.68 

1.73 

14.6 

I03A 

See Jackson 2011 

5 

6 

3 


Table 1. Basic information for the systems studied, including the redshifts of lens and source (where known), maximum separation of 
the lensed images, observation time and frequency, and achieved noise level in the maps. References are given to the discovery papers for 
each lens, and to the sources for the measurements of the redshifts. In each case the on-source integration time is approximately 75% of 
the total exposure time. Key to references: W02 = Wisotzki et al. 2002, M05 = Morgan et al. 2005, R02 = Reimers et al. 2002, B97 = 
Bade et al. 1997, B98 = Burud et al. 1998, KOO = Kneib, Cohen & Hjorth 2000, I03A = Inada et al. 2003a, I03B = Inada et al. 2003b, 
E06 = Eigenbrod et al. 2006, 006 = Ofek et al. 2006. 


Object 

Phase calibrator 

Date 

Exposure 

Frequency 

HS 0810+2554 

JVAS 0813+2435 

31/03/2014 

8h 

1287-1799MHZ 

RXJ0911+0551 

SDSS J0901+0448 

26/04/2014 

8h 

1287-1799MHZ 


Table 2. Details of the e-MERLIN observations of two of the sources. 


package. The data were fringe fitted to remove instru¬ 
mental delays using the phase calibrator observations, 
and a phase and amplitude solution was performed us¬ 
ing the phase calibrator and a point source model. Maps 
were also made of the phase calibrator and used where 
necessary to iterate the phase and amplitude calibra¬ 
tion, and the flux density solution was transferred from 
the flux calibrators. The resulting calibration was then 
applied to the target sources, which were imaged using 
natural weighting in order to achieve the best possible 
signal-to-noise. In general we obtain image noise levels 
close to the theoretical value, usually about 3/iJy/beam 
for 6 hours of on-source time. 


2.2.2 e-MERLIN observations 

Two of the objects (HS 0810+2554 and 
RX J0911+0551) were also observed with the e- 
MERLIN array (Table 2). The observations were 
carried out at L-band, with a bandwidth covering 
the wavelength range 1287-1799 MHz. In addition 
to the target sources, observations of nearby phase 
calibrators were carried out, with a cycle of 7 minutes 
on the target and 3 on the phase calibrator. Additional 
observations of 3C286 were carried out in order to set 
the flux scale, and of the bright point source OQ208 
in order to calibrate the bandpass. Data reduction 
followed standard procedures (Argo 2015) including 
fringe fitting to all calibrator sources to determine 
delays, phase and amplitude calibration using the 
nearby phase calibrator, and determination of the 
flux scale and bandpass calibration, with allowance 
for the spectral index of the calibrator. The telescope 
weights were modified using standard values for L-band 
provided by the observatory, and final maps were made 
in the AIPS software package distributed by NRAO. 
Mapping in the case of HS 0810+2554 was complicated 
by the presence of a 200-mJy confusing source 6' from 
the target; this source was mapped simultaneously 
with the target, and was also used to refine the phase 
calibration of these observations. Noise levels achieved 


Object 

Type 

Gpt. 

Flux density 
(radio,/iJy) 

HE 0435-1223 

Gross 

A 

36.0+2.1 



B 

26.4+2.1 



G 

34.3+2.1 



D 

16.1+2.1 

HS 0810+2554 

Fold 

A 

85.1+3.7 



B 

83.7+3.7 



G 

60.0+3.7 



D 

49.1+3.7 

RX J0911+0551 

Gusp 

A 

26.9+2.2 



B 

53.2+2.2 



G 

19.7+2.2 



D 

9.4+3.0 



G 

18.3+2.2 


Table 3. Radio flux measurements for the sample of four-image 
lens systems observed with the VLA/e-MERLIN, in which fluxes 
can be measured. For HS 0810+2554 the flux densities are from 
the VLA map at 8.4 GHz; the corresponding flux densities in the 
e-MERLIN image are (161,173,129,216) for A,B,C and D respec¬ 
tively, with errors of approximately 30p,Jy in each case. 

in these observations were about 15-30/iJy, depending 
on the details of the mapping strategy. 


3 RESULTS AND MODELS 

All four objects were detected in these observations, 
of which all but SDSS J0924+0219 have individually 
measured flux densities for each lensed image. We dis¬ 
cuss the results for each object separately, before mak¬ 
ing more general remarks about the measurements. The 
radio flux densities are given in Table 

The approach to modelling the observations is the 
same in each case. First, we make a preliminary assess¬ 
ment of whether the radio map is consistent with lensing 
of a point source into point images. We do this by mod¬ 
elling the structure, in each case, with four point-spread 
functions (PSFs; Table [3), whose extent is known accu¬ 
rately from the radio CLEAN procedure. In this model, 
we fix the separation between the four individual com¬ 
ponents using measurements from archive HST images 
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as reported by the CASTLES astrometry (which are 
accurate to a few milliarcseconds), but the overall reg¬ 
istration of the image has been allowed to vary. There 
are thus six free parameters in the model, two for the 
registration, and four from the flux densities of the indi¬ 
vidual points. Second, we make a lens model using con¬ 
straints from the image-plane radio map; for this, we 
assume a singular isothermal mass distribution (except 
in the case of HE 0435—1223) for the lenses together 
with a contribution from external shear. The source is 
assumed to be of Gaussian profile, and the resulting 
image plane is compared to the data, optimising the 
lens galaxy parameters together with the source posi¬ 
tion, flux density, size and ellipticity. We note that mod¬ 
ern wide-bandwidth interferometers at centimetre wave¬ 
lengths, such as the VLA and e-MERLIN, come close 
to filling the u — v plane. Because of this virtually filled 
aperture, there is no need to fit the data directly in the 
u—v plane. This contrasts with the situation in early At¬ 
acama Large Millimetre Array (ALMA) datasets used 
to map sub-millimetre lenses (e.g. Hezaveh et al. 2013). 

3.1 HS 0810-h2554 

HS 0810-1-2554 was discovered by Reimers et al. (2002) 
and consists of four images with the two southwestern, 
bright images in a merging pair configuration. In HST 
imaging (Reimers et al. 2002) the leasing galaxy is de¬ 
tected, with an unknown redshift (it is estimated as 0.89 
by Mosquera & Kochanek 2011 from the separation and 
the redshift distribution of existing lenses). These im¬ 
ages also show a 0.7-magnitude difference in brightness 
between the components of the merging pair, contrary 
to the expectations of simple models, but this is likely 
to be due to microlensing in the lens galaxy. The source 
is a narrow absorption line quasar, with relativistic out¬ 
flows detected using X-ray absorption spectra (Chartas 
et al. 2014). These high velocity outflows may be mag¬ 
netically driven. 

Our radio maps from the VLA and e-MERLIN are 
shown in Eig. The components in the VLA image 
appear extended, and a faint arc is visible around the 
bright merging pair. This extension can be quantified 
by attempting to model the lensed structure only with 
point sources, with separations fixed to that of archival 
HST data. The best fit shows significant residuals, in 
particular around the bright merging components, but 
also at a lower level around the line connecting images 
B and C (Eig. 1^. It is obvious from visual inspection 
that the shape of the A-B complex in the data is signifi¬ 
cantly different from that of a two-Gaussian realisation. 
We therefore conclude that the source is likely to be ex¬ 
tended and model it as such. 

In order to model the extended source, we have as¬ 
sumed a simple lens model (singular isothermal ellipsoid 
plus external shear) and treated the source as an ellipse 
with a flux density, axial ratio, position and orientation 
which are free to vary. Eor each iteration of the model, 
the source is projected through the lens model, and the 
result convolved with the PSE of the radio map. Min¬ 
imisation of the between the model image and data 
is carried out, where the x^ value is determined from 


a comparison of the model with the image in all re¬ 
gions where either model or image contains flux above 
2(7. Gorrelations between pixels are neglected. The po¬ 
sition of the galaxy is fixed by the use of the Hubble 
Space Telescope (HST) image during this process. In 
practice, the quality of the fit does not depend signifi¬ 
cantly on this quantity, provided that the source is al¬ 
lowed to move to keep the same distance between it 
and the galaxy. An acceptable fit of x^=1-6 is obtained 
with such a procedure; the parameters of this fit are 
given in Table and shown in Eig. and ^ The two 
close images, A and B, are of approximately the same 
flux density, as expected in the absence of millilensing- 
scale substructure and in contrast to the measurements 
in the optical and near infrared. The model implies a 
magnification of about 25 for the brightest image, yield¬ 
ing an intrinsic flux density of 3.5 /iJy for the source. 
The implied magnification is a factor of 2 less than that 
of the model by Assef et al. (2011), but HS 0810+2554 
is in the high-magnification regime where the source is 
contained within a very small astroid caustic, and mi¬ 
nor movements in the source position produce major 
changes in the implied magnification. 

In order to derive uncertainties, the parameters 
have been run through a Markov Chain Monte Carlo 
(MCMC) process using the publicly available emcee 
routine (Eoreman-Mackey et al. 2013). We have as¬ 
sumed a number of hard priors, namely limits of 0"- 
1" for the Einstein radius of the galaxy, limits of 0 < 
7 < 0.5 for the external shear contribution, a require¬ 
ment that the source flux and source width are positive, 
and that ellipticities of the source and galaxy must be 
<1. Einally, we have imposed a Gaussian prior on the 
position of the lens galaxy, based on the measured po¬ 
sition in the HST image and with a width of 10 mas. 
The results (Table make it clear that the source is 
extended by approximately 10 mas in the source plane, 
corresponding to about 70 pc in physical scale, at a po¬ 
sition angle of approximately 50°. As usual with strong 
lens systems, we obtain a relatively accurate measure¬ 
ment of the Einstein radius of the lens galaxy, which 
is modelled as being almost circular. This is consistent 
with its appearance on archival HST images. 

Although the existing e-MERLIN images do not 
have very high signal-to-noise, they do allow us to mea¬ 
sure an approximate overall spectral index, because the 
resolution of the e-MERLIN 1.6-GHz image is very sim¬ 
ilar to that of the VLA at 8.4 GHz. This spectral index 
is moderately steep, at —0.55 ±0.1. 

3.2 HE 0435-1223 

HE 0435—1223 was discovered by Wisotzki et al. (2002) 
and identified as a four-image system with an early-type 
lens galaxy. The lens redshift was obtained by Morgan 
et al. (2005), who also found that the lens is part of a 
group of galaxies. Microlensing was detected in a sub¬ 
sequent monitoring campaign (Kochanek et al. 2006) 
which probably affects the A component most strongly 
(Ricci et al. 2011, Courbin et al. 2011), and it has also 
been shown that the broad-line region in this object 
is subject to microlensing (Sluse et al. 2012; Braibant 
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Figure 1. Left: VLA greyscale radio map of HS 0810+2554 at 8.4 GHz. The beam is of full width at half maximum (FWHM) 300x240 mas 
in position angle —65.17°. Archival HST Near-Infrared Camera and Multi-Object Spectrometer (NICMOS, red) and Advanced Camera 
for Surveys (ACS, green) contours have been performed using image C. The conventional nomenclature of the images (Reimers et al. 
2002) is that the merging complex in the southwest consists of images A and B, with A being further north. Right: e-MERLIN image of 
HS 0810+2554 at approximately the same resolution, but a frequency of 1.6 GHz. The noise level is approximately 29)nJy/beam; all the 
images are detected at > 5a significance. 
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Quantity 


HS0810+2554 HE0435-1223 RXJ0911+0551 


Source position/mas 
Source FWHM along major axis/mas 
Source flux/p,Jy 
Source b/a 

Source position angle 
Galaxy critical radius/mas 
Density slope (2 = isothermal) 
Galaxy ellipticity 
External shear 


0.1 E, 13.0 S 
12+1 
3.6+0.2 

(47±5)° 

473+10 

= 2.0 

0.0003+0.0003 
0.023+0.006, (29+4)° 


-68 E, 18 S 

sotl 

2 9+0.3 


^-0.4 

= 1.0 


113816 ^ 

z.uu_o 06 

0.26+0.02 

0.039l°;°°i, (-30±7)° 


468 E, 28 S 

3 7+0.3 


' - 0.2 
= 1.0 


1047Z 


7l38 
=2.0 
n 1 cr+O.OO 
*^•^^-0.09 

0.373l°°fi, (9±2)° 


Table 4. Model fitting results for the three lens systems for which lens modelling is possible (all observed lenses except SDSS 0924+0219). 
The source position is quoted relative to the galaxy position, and all distances are given in units of milliarcseconds. For HE0435—1223 
and RXJ0911+0551, the galaxy critical radius corresponds to the Einstein radius measured along the major axis. 



Figure 2. Model of HS0810+2554, using a point-source model fit (see text). The data, model, and residual are plotted; unlike the case 
with the extended source model, there appear to be significant features in the residual. Note that in this case, and unlike the case of 
the extended source model, the overall shape of the A-B complex is not well reproduced. The bar in each panel represents 1", and the 
colour scales in the sidebars are in units of p,Jy/beam. 
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Figure 3. Models of HS0810+2554, using a singular isothermal sphere model plus external shear (see text) together with a best-fit 
extended source. The reconstructed source, data, model, and residual are plotted. The white bar represents 10 mas in the panel of the 
reconstructed source, and 1" in all other cases. In these and subsequent figures, the data and model are plotted on the same colour scale, 
and the residuals are scaled to the minimum/maximum of the residual map. The colour-bars represent flux densities in p,Jy/beam; the 
colour scale of the source is arbitrary, but its parameters are given in Table A good fit is obtained with an unlensed source size of 
between 10-15 mas. 



Figure 4. MCMC realisations of the model of HS 0810+2554, 
plotted as probability density of source size against source axis 
ratio. The preferred source size corresponds to a linear scale of 
approximately 100 pc, a unique resolution for such a faint radio 
source. 


et al. 2014). Modelling of the lens is able to reproduce 
well the positions of the lensed images, using only a 
single-galaxy deflector model (Sluse et al. 2012). Fadely 
& Keeton (2012) examined and modelled this object ex¬ 
tensively in a search for evidence of substructure in the 
lensing galaxy, using their L'-band mid-infrared fluxes 
of the four components as inputs for the modelling. 
Since much of the quasar mid-infrared emission comes 
from a circumnuclear torus, this may be immune to mi- 
crolensing as the torus is likely to be relatively large. 

Our radio map (Fig. shows all four lensed im¬ 
ages, at positions negligibly different from those of the 
optical and infrared HST images (obtained from the 
HST archive). Once again, therefore, we have modelled 
the radio map assuming that it consists of four point 
sources, whose separation is determined by the HST 


Figure 5. VLA radio map of HE 0435—1223. The greyscale runs 
from 0 to 20p,Jy/beam. The beam is circular, and of FWHM O'.Mh. 
Archival NICMOS (red contours) and ACS (green contours) im¬ 
ages are overlaid. Registration of these images has been performed 
using image A. The conventional nomenclature of the images is 
that A is the easternmost and B,C,D proceed clockwise around 
the lens galaxy. 


optical image. The results of this procedure are shown 
in Fig. There are hints of emission outside the four 
point sources, particularly close to image C, but these 
are at the 2-cr level and would need deeper observa¬ 
tions to confirm or rule out. The image fluxes are given 
in Table and, with a ratio of 1.05:0.77:1.00:0.47 be¬ 
tween A:B:C:D images, differ significantly from the L' 
ratios 1.71:0.99:1.00:0.81 of Fadely & Keeton (2012). In 
particular, the A/C and C/D ratios differ by about 3cr 
between the two sets of observations, the difference in 
C/D ratio being particularly obvious from Fig. 

This difference in flux density ratios, and its ex¬ 
planation, is a difficult and intriguing problem. Fadely 
& Keeton (2012) undertook extensive modelling of this 
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Figure 6. Models of HE 0435—1223, using a point-source model fit (see text). The data, model, and residual are plotted; unlike the 
case with the extended source model, there appear to be significant features in the residual. The bar in each panel represents 1". The 
numbers on each colourbar are in units of /rJy/beam. 
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Figure 7. Models of HE 0435—1223/ From left to right: source model; data; image-plane model; residual. The white bars represent 
100 mas in the source plane (left panel) and 1 arcsecond in all other panels. Three different fits are shown: (top) small 3 mas source, 
(middle) optimal 80 mas source, (bottom) large 200 mas source. Although the 80 mas source is preferred by our MCMC analysis, smaller 
sources provide an equally good visual fit, while larger sources lead to significant model residuals. The colour-bars represent flux densities 
in /rJy/beam; the colour scale of the source is arbitrary, but its parameters are given in Table ^ 




n 



system using a softened power-law for the primary 
galaxy, together with a singular isothermal model for 
the nearby galaxy, G22, which is about four arcsec- 
onds to the SW. They also included a shear compo¬ 
nent, to account for the more general shear field of the 
cluster. Smooth models with a slightly shallower slope 
than isothermal failed to reproduce the infrared fluxes, 
and further analysis showed that the Bayesian evidence 
favoured a contribution due to substructure around A. 
The observation driving this conclusion was the A/C 
ratio, which smooth models preferred to be between 1.4 
and 1.5, as opposed to the higher value in the infrared 


data. Our radio data, on the other hand, prefer a much 
lower value for the A/C ratio, together with a much 
fainter D component. Since the radio is almost certainly 
not affected by microlensing, this is a puzzling result. 

In our next step of modelling the data, we assume 
the source is point-like and include the observed time 
delays (Courbin et al. 2011) as additional constraints. 
Modelling the lens as an ellipsoidal power-law with ex¬ 
ternal shear and a SIS at the position of C22 (z = 0.78, 
Chen et al. 2014), realistic, smooth models are unable 
to reproduce the data. The best model, further con¬ 
strained using strong Caussian priors on the positions 
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(lcr =10 mas) of the galaxies, agrees with those observed 
time delays within the errors but not with the image 
fluxes, yielding a of 14.9 for four degrees of freedom 
In particular, the predicted B/C and D/C ratios are in 
disagreement at the 1 . 3 cr and 3 cr levels, respectively. For 
the main lensing galaxy, the best model prefers an el- 
lipticity of e ~ 0.28 and power-law slope of 7 ' ~ 2.24. 
We note that the steep density slope may be driven by 
the model trying to fit the flux ratios. 

The flux ratio anomaly seen with the smooth model 
could be explained by invoking substructure. To explore 
this possibility, we follow the approach of Fadely & Kee¬ 
ton ( 2012 ). Briefly, substructure clumps are modelled 
using a pseudo-Jaffe profile, and a wide range of masses 
are considered, whose masses enclosed within the Ein¬ 
stein radii range from 10^ — 10^ Mq. Modelling the 
smooth lens component as before, we find that clumps 
placed near images A, B, or C do not improve the fit. 
However, clumps over a large range of masses placed 
near image D can bring the model into good agreement 
with the data, yielding a of 1.5 for one degree of 
freedom.^ Like Fadely & Keeton ( 2012 ), we find that 
more massive clumps can be placed farther away from 
the image, while less massive clumps must be placed 
nearer. Furthermore, the clumps cannot he within ap¬ 
proximately two Einstein radii of image D. Otherwise, 
image splitting can magnify the image, making matters 
worse. As for the main lensing galaxy, ellipticities of 
e ~ 0.33 and steeper-than-isothermal power-law slopes 
of 7 ' « 2.33 are preferred. 

Alternatively, if we do not invoke substructure, an¬ 
other possible solution arises if the radio emission region 
is extended and differentially magnified (see e.g. Ser¬ 
jeant 2013). Because the size of the caustic is less than 
approximately 400 mas, a wide range of source sizes be¬ 
low this scale can reproduce the data. We again use a 
two-deflector model, including the main lensing galaxy 
as an ellipsoidal power-law with a contribution from ex¬ 
ternal shear and G 22 as a SIS. The source is modelled 
as a spherical Gaussian. Owing to the large number of 
image pixels (3150 pixels) used to constrain the model, 
including time delay constraints for a point source at 
the position of the source does not significantly affect 
our results, and so we include them for consistency with 
the previous analyses. Additionally, to try to account for 
the noise correlation in the data and to be conservative 
in our parameter inference, we follow the suggestion of 
Riechers et. al. (2008) and increase the input noise level 
by a factor dependent on the noise correlation length 
scale (a factor of three, in this case). Our best model 
achieves a y^ of 854 for 3141 degrees of freedom.^ Af¬ 
ter marginalizing over all lens model parameters, our 
MGMG analysis finds that an isothermal slope is pre¬ 
ferred for the main lensing galaxy ( 7 ' = 2 . 00 ^o oe) 


^ Because of the intrinsic nonlinearity of the lens modelling and 
the strong Gaussian priors placed on the galaxy positions, calcu¬ 
lating the the number of degrees of freedom is nontrivial. Thus, 
we have chosen to count each galaxy position parameter as half a 
degree of freedom, and so the “true” number of degrees of freedom 
may be ±2 of the number we quote here. 


that the source is of EWHM 8 OI 5 mas. This optimal 
source size leads to an image-plane model that shows 
discernible extended structure (Eig. [^. By visual in¬ 
spection, we find that source sizes more than an order 
of magnitude smalleij^ can also fit the data reasonably 
well but result in point-like images at these resolutions 
and leave larger model residuals. Sources a factor of 
two larger, on the other hand, are clearly unfavorable 
by both visual inspection and the MGMG analysis. 

Of the possible choices for explaining the data, we 
prefer the option that the source is extended and dif¬ 
ferentially magnified. As radio sources are likely to be 
more extended than their optical or mid-infrared coun¬ 
terparts, this seems to be the most natural choice. Eur- 
thermore, finite source size effects would likely be re¬ 
quired in order to simultaneously explain the flux ratio 
anomalies in the mid-IR as well. We note, however, that 
due to the non-Gaussianity and correlation of the noise 
in the immediate regions surrounding the lensed images, 
we do not strictly limit our conclusions by the results of 
the MGMG analysis. Instead, we provide a conservative 
upper limit of 200 mas for the source size. Eor sources 
larger than this, a visual inspection of the model resid¬ 
uals clearly shows that the source has become too large. 

3.3 RX J0911+0551 

RX J0911+0551 (Bade et al. 1997) is a cusp- 
configuration lens system, with three close images (A,B 
and G) and a fourth image some distance to the west. 
The mass environment is relatively complicated; the 
lens lies close to a massive cluster about 40" away and 
in addition to the primary galaxy, a second galaxy lies 
close to the system, complicating the process of mass 
modelling. Our VLA 5-GHz image is shown in Eig. 
and has a noise level of just under 2 /iJy/beam. All four 
lensed images are clearly detected in the radio map, and 
in addition we detect radio emission at the position of 
the lensing galaxy. We do not detect any of the com¬ 
ponents in the e-MERLIN 1.5 GHz image, which has a 
noise level of 16/iJy/beam. 

To quantify the non-detection, and thus derive lim¬ 
its on spectral index between the e-MERLIN and VLA 
observations, the e-MERLIN observation was conserva¬ 
tively simulated with four components of the size of the 
VLA beam (^500 mas), whose flux densities were given 
by the VLA map. Gaussian noise was added to the map 
at the same level as the observations (i.e. RMS 16/iJy). 
The components were used to generate a u — v dataset 
with the sampling function and noise level of the ac¬ 
tual e-MERLIN observations. This was imaged and the 
lower limit on the spectral index resulted from the in¬ 
put fluxes for which the components could just not be 
recovered from the simulated image by model fitting. 
The lower limit on the spectral index was found to be 
a = —0.5. 

We have again attempted to model the image plane, 
represented by the VLA map, as a sum of point spread 

The image plane is appropriately oversampled to ensure that 
fluxes are calculated accurately. 
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Figure 8. JVLA radio map of RX J0911+0551 (greyscaled from 
0 to 20p,Jy) with contours from archival HST/NICMOS data 
superimposed, aligned on image D. The beam is of FWHM 
560x390 mas in position angle —39.2°. The three close images 
to the left are A,B,C (from south to north), and image D is at 
the right of the picture. The lensing galaxy (between the images) 
is radio-loud, with a flux density of about 16/rJy. 


Figure 10. VLA radio map of SDSS J0924+0219, with greyscale 
from 0 to 20)11 Jy/beam, and a restoring beam of 696x416 mas 
in PA —43°. The object is detected, but individual flux densities 
for the images are impossible to measure. Image registration to 
archival HST/NICMOS data (contours) has been done by eye, 
but this procedure is not well-deflned to better than the absolute 
astrometry of the HST image. 


functions whose relative position is determined by the 
HST astrometry (Fig.|^. Here it does appear that point 
models have difficulty in reproducing the structure, in 
particular the shape of the A-B-C complex, although 
we recover good estimates for the flux densities of each 
image (TablejS). The flux densities of the A,B and C im¬ 
ages in the r^io have a ratio very close to 1:2:1, close 
to that expected by the cusp relation and suggesting 
that substructure does not need to be invoked in this 
case. This contrasts with the optical flux densities, and 
in particular with the flux density ratio 2:2:1 between 
A, B and C measured by Burud et al. (1998). Hence the 
optical measurements are almost certainly affected by 
microlensing. By contrast, Sluse et al. (2012) model this 
lens using astrometric constraints and a model consist¬ 
ing of a singular isothermal ellipsoid plus shear. They 
obtain image flux ratios (A:B:C:D) of 1:1.87:0.88:0.34. 
Our corresponding values are 1:2.05:0.73:0.35, consis¬ 
tent within the errors with Sluse et al.’s model. The 
source flux density predicted by this lens model is about 
5/iJy. 

Motivated by the residuals observed in the point- 
source model, we have again fitted a model in which 
a Gaussian-shaped extended source is lensed. Provided 
that the source size is not very small, neither it nor the 
source shape is well constrained (Table |^. The mod¬ 
elled shear is large, suggesting that we are seeing the 
influence of the cluster to the south. The good overall 
fit to the data, x^=3518 for 2290 degrees of freedom, 
gives no significant evidence for effects of substructure 
in the lens galaxy. 


3.4 SDSS J0924+0219 

SDSS J0924+0219 (Inada et al. 2003b) is a very inter¬ 
esting object because it has such a large optical flux 
anomaly, almost certainly as a result of microlensing 
(Morgan et al. 2008). This anomaly decreases with in¬ 
creasing observation wavelength (Floyd, Bate & Web¬ 
ster 2009) allowing the size of the accretion disk to be 
measured, but the anomaly persists over a period of 
^7 years (Fame et al. 2011) raising the possibility that 
some of it may be due to the effects of lensing by sub¬ 
structure. 

Although we detect the object at a reasonable level 
of significance (Fig. [^, we unfortunately do not have 
sufficient signal-to-noise, in the three hours of obser¬ 
vation time allocated to this object, to measure the 
flux densities separately. At a total flux density of 
^ 15)Ujy, this is by far the weakest of the objects stud¬ 
ied. Using an isothermal model, plus external shear, fit¬ 
ted to the positions reported by Inada et al. (2003b), 
we obtain magnifications for the four components (A, 
B, C and D) of approximately 13, 5, 5 and 11. Our 
overall flux density of ^ 50 )Ujy implies an unlensed 
source flux density of about 1.5)Ujy. This is the sec¬ 
ond faintest radio source yet detected, the faintest being 
SDSS J1004+4112 (Jackson 2011); further observations 
of the sample of radio-quiet quasars are likely to yield 
the first detected nano-Jy radio source. 
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Figure 9. Models of RXJ 0911+0551, using a point-source model fit (above, see text). The data, model, and residual are plotted. In 
this case, the point-image model appears to have significant difficulty in fitting the A-B-C complex at the eastern end of the system. The 
bar in each panel represents 1". Model using an extended source (below). The panels contain the source plane, the data, the model and 
the residual. The area around the galaxy has been blanked and excluded from the fit. In all cases the numbers on the colourbars are in 
/rJy/beam; the source panel colourscale is arbitrary, but the source parameters are given in Table 


4 DISCUSSION AND CONCLUSIONS 

4.1 Radio properties of the lensing galaxies 

We detect the lensing galaxy in only one of the four 
objects: RX J0911+0551, at a level of 18/iJy. The radio 
flux density of the lensing galaxies in the other three 
objects is < 6 — 8 /iJy (3cr). 

The lensing galaxy flux density in RX J0911+0551 
corresponds to a luminosity of 5 x 10^2 WHz-\ which 
is at the top of the range that can plausibly be at¬ 
tributed to star formation; this range is bounded by 
the 10^1 WHz-i of the Milky Way and lO^^ WHz-^ 
for extreme star-forming galaxies. It is also close to 
the boundary between star-forming radio emitters and 
AGNs found by Kimball et ah (2011) in their study of 
radio luminosity functions of nearby {z ^ 0 . 2 ) AGN. 
There is no evidence for ongoing star formation in the 
lens galaxy, which seems to be an early-type galaxy, in 
other wavebands. For example, Burud et al. (1998) con¬ 
ducted optical and near-IR imaging and found that the 
lensing galaxy has a similar red colour to that of nearby 
members of the cluster of which it is part. 


4.2 The nature of radio-quiet quasars 

This work has resulted in the measurement of radio flux 
densities, and in some cases resolved radio structure, in 
a significant number of radio- “quiet” lensed quasar sys¬ 
tems. Such radio sources, if unlensed, would be beyond 
the reach of current instruments in all but exceptional 
observing times, and are objects whose study will only 


become routine with the SKA. The luminosity of an ob¬ 
ject of intrinsic flux density of 1/iJy and flat spectral in¬ 
dex is about l.OxlO^^WHz”^ at ^ = 0.5, 5xlO^^WHz“^ 
qX z = 1 and 30xl0^^WHz“^ at 2 : = 2 , orders of mag¬ 
nitude below what is typically accessible with current 
surveys except at low redshift (e.g. fig. 4 of Condon et 
al. 2013). Previous studies of radio-quiet quasars have 
focused on optically bright quasars, such as the Palomar 
Green sample (Kellermann et al. 1994). In accordance 
with the radio-optical correlation noticed by White et 
al. (2007), these objects have typical radio flux densities 
of a few hundred /iJy, two orders of magnitude brighter 
than the intrinsic flux densities of the objects studied 
here. 

We have observed four of the 15 known optically- 
selected, four-image quasar lenses with S > — 20 °, and 
all of them have intrinsic radio flux density of between 
1 and 5/iJy. Of the other nine, three are known to 
have significant radio emission. PGl 115+080 has VLA 
archival data at 8.4 GHz taken in the compact (D) 
configuration, which yield a total radio flux density of 
153T17 /iJy, although the resolution of a few arcseconds 
does not allow the flux density of individual components 
to be determined. However, the likely magnifications in 
this lens system suggest that the intrinsic flux density 
of the source is also a few /iJy. A similar result can 
be derived for the lens system RX J1131 —1321, which 
was found to have significant radio emission by Wuck- 
nitz & Volino (2008). Finally, H1413+117 is a radio¬ 
intermediate object which has been studied with the 
VLA by Kayser et al. (1990). Further lenses from the 
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Figure 11. Intrinsic radio flux densities versus intrinsic I-band 
magnitudes for a sample of optically-selected quasar lenses. Both 
radio and optical flux densities have been demagnifled using lens 
models. Data and models have been derived from Ratnatunga, 
Griffiths & Ostrander 1999, Wisotzki et al. 2002, Reimers et al. 
2002, Burud et al. 1998, Inada et al. 2003 a,b, Ghosh & Narasimha 
2009, Anguita et al. 2009, Jackson 2011, Wucknitz & Volino 2008, 
Kayo et al. 2007, Oguri et al. 2008, Assef et al. 2011, in some case 
supplemented by further modelling. The locus of radio vs. optical 
flux densities reported by White et al. (2007) is sketched near the 
top right. 


COSMOS survey (Faure et al. 2008, Jackson 2008) do 
not have significant radio emission (Schinnerer et al. 
2007, Schinnerer et al. 2010) in the VLA-COSMOS sur¬ 
vey, and three quad lens systems from the Sloan Quasar 
Lens Search (SDSS J1138+0314, SDSS J1251+2935, 
SDSS J1330+1810) do not yet have deep radio imag¬ 
ing. It is therefore likely that at least half of optically- 
selected quasar lens systems will show radio emission 
at the micro-Jansky level, if examined carefully, and a 
more complete census will be the focus of future work. 

The existing data are plotted in Fig. El The radio 
fluxes have been derived from the literature (see the fig¬ 
ure caption), with in some cases a limit of 1 mJy inferred 
from their absence from the FIRST 1.4-GHz catalogue. 
The current sample is small. However, we note that 
the median radio flux density inferred by White et al. 
(2007) from their stacking analysis is about 50-70 /iJy 
at 20 < / < 21. This, combined with the distribution 
of detections in our optically-selected sample, suggests 
that there is a large scatter in radio flux densities at 
this optical magnitude, if not an outright bimodality. 

There are a number of theoretical models for the 
origin of radio emission in radio quiet quasars, each 
of which makes rather different predictions for what 
should be observed. The first possibility is that of a 
smaller version of a radio-loud quasar, where a flat- 
spectrum radio core and steeper-spectrum radio jet 
emission are present in some ratio (Urry & Padovani 
1995, Ulvestad et al. 2005). In this case, we would ex¬ 
pect steep-spectrum emission on scales of parsecs to 
tens of parsecs in addition to a compact, flat-spectrum 
radio core. A second possibility is the emission from ra¬ 
dio starbursts in a similar manner to nearby examples 


such as M82 (e.g. Condon et al. 2013). This emission 
is expected to be optically thin synchrotron from su¬ 
pernova remnants, but should extend over most of the 
galaxy disk and therefore have a characteristic size of 
about 1 arcsecond, or 5-10 kpc; studies of star-forming 
radio sources in the HDF, with mean redshifts around 
1, have shown that the radio emission nearly always dis¬ 
plays a similar characteristic size (Muxlow et al. 2005). 
Alternatively, the radio emission could be produced by 
thermal processes close to the accretion disk. Sugges¬ 
tions for this include optically thin free-free emission 
from a disk wind (Blundell & Kuncic 2007) and emis¬ 
sion from magnetically-heated coronae (Laor & Behar 
2008). In both cases the emission would be expected to 
originate close to the centre. In the case of disk winds, 
this scale is likely to be at least 0.1-lpc, but for magnet¬ 
ically heated coronae the scale would be smaller. In both 
these cases, however, the radio spectral index should be 
approximately flat. 

We can use our data to confront the models in a 
number of ways. The first relevant result is the source 
sizes implied by our lens modelling, which are of the or¬ 
der of 70 pc in HS 0810+2554. In HE 0435-1223 we find 
a characteristic size which is more uncertain, but prob¬ 
ably greater than a few milliarcseconds and certainly 
less than 200 mas (<2 kpc). In RX J0911+0551 we 
again detect a significantly extended source, but whose 
size is likely to be 100-150mas (^1 kpc) rather than 
the larger values which would be expected for a typical 
star-forming region in a radio source. Either we have an 
unrepresentative sample, or the star-formation model is 
disfavoured compared to the non-thermal AGN hypoth¬ 
esis. This is in contrast to the inference of star-formation 
as the cause of the radio emission, made by Wucknitz & 
Volino in the case of RX J1131 —1231. A definitive test 
should be available using VLBI, as a non-thermal source 
should contain high-brightness emission at the /iJy level 
which is detectable with current VLBI sensitivities. 

A second, although rather more equivocal, result 
concerns the measured spectral indices. The spectral in¬ 
dex of the radio emission in HS 0810+2554 appears to 
be steep, consistent with synchrotron emission from ei¬ 
ther a non-thermal source or a star-forming component. 
Both the spectral index and the characteristic size dis¬ 
favour coronal emission and emission from disk winds, 
which would be expected to be relatively flat-spectrum 
and to be emitted from a smaller region. However, the 
VLA and e-MERLIN observations of RX J0911+0551 
may be consistent with such models. Because the spec¬ 
tral index limit is currently relatively loose {a > —0.5), 
further e-MERLIN observations are needed to make a 
more definite statement. 

4.3 Substructure in lensing galaxies 

Radio flux densities in four-image lens systems are im¬ 
portant because they potentially give indications of sub¬ 
structure in lensing galaxies (or along the line of sight), 
in the form of flux anomalies (Mao & Schneider 1998; 
Dalai & Kochanek 2002). In particular, violations of 
the cusp and fold relation allow us to quantify the lev¬ 
els of substructure present (e.g. Xu et al. 2009, Xu et 
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al. 2015). Flux ratios between images can also be af¬ 
fected by propagation effects (scattering in the case of 
the radio waveband), microlensing (in the case when the 
source is smaller than the intrinsic size of the Einstein 
radius of stars in the lensing galaxy, or about 1 /rarcsec), 
variability (in the case of a source which varies signif¬ 
icantly over the time delay between the images) and 
source size (which can affect flux ratios, e.g. Amara et 
al. 2006; Metcalf & Amara 2012). The use of cusp and 
fold relations, rather than the observation of disagree¬ 
ment with a smooth model, is important because other¬ 
wise the effect of substructure on the image flux ratios 
can be partially absorbed by varying the smooth model. 

We detect no new evidence in our objects for sub¬ 
structure. Indeed, the flux ratios of RX J0911+0551 
obey the cusp relation within the errors of the mea¬ 
sured radio flux densities. The flux ratios of the fold sys¬ 
tem HS 0810+2554 are also consistent with a smooth 
model, as the brightnesses of the merging images are 
equal to within the errors, unlike the optical case in 
which microlensing affects one of the images. The case of 
HE 0435—1223 is more interesting. Here we obtain flux 
density ratios which are inconsistent with not only the 
optical, but also the mid-IR bands. There are a number 
of explanations for this which we do not favour. Radio 
microlensing could affect the radio flux densities, but 
would require a very small radio source size (microarc- 
seconds rather than milliarcseconds). Variability of the 
radio source, together with a time-delay, is also possi¬ 
ble. However, intrinsic variations in typical radio-faint 
quasars are on timescales of several weeks to months 
with fluctuations at the 10 — 20% level (Barvainis et 
al. 2005). In order to reproduce these observations, a 
variation in the source flux of ^40% would be required 
within the time-delay scales measured by Courbin et al. 
(2011) (—6.5 and —14.3 days for B-D and C-D respec¬ 
tively). Elux anomalies due to variations in the mid- 
IR are still less likely, because the L'-band is expected 
to have a significant contribution from a dusty torus, 
whose size is > Ipc, implying a light crossing-time of 
several years or more. Differential extinction due to dust 
at non-radio wavelengths is not a likely culprit, as the 
colours (Wisotzki et al. 2002) and the continuum slopes 
(Morgan et al. 2005; Wisotzki et al. 2004) are nearly 
identical for the four lensed images. While substructure 
can be used to explain the various flux ratio anoma¬ 
lies seen in the optical, mid-IR, and radio, it may prove 
difficult and would perhaps require fine-tuning to simul¬ 
taneously explain all the observations with substructure 
alone. Instead, at least some of the explanation is likely 
to be the effects of finite source sizes; especially given 
the results in the other objects, an intrinsic radio source 
size of order a few parsecs is the explanation that we 
favour. 
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